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O ' Abstract 
00 ' 

We study the mechanism leading to the metallization of the /3-SiC(001) Si-rich surface induced 
by hydrogen adsorption. We analyze the effects of band bending and demonstrate the existence of 

X 

a quasi-2D electron gas, which originates from the donation of electrons from adsorbed hydrogen 
to bulk conduction states. We also provide a simple model that captures the main features of the 
results of first-principles calculations, and uncovers the basic physics of the process. 
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FIG. 1: Band structure of the model of the metallized surface (shown schematically in the inset). 
The broken line is the Fermi level. Circles at T and J indicate the quasi-2D bands. 

The Si-rich /3-SiC(001) 3x2 surface has been reported to become metallic when exposed to 
atomic hydrogen |2J. This is quite surprising, as H-saturation of semiconductor surfaces 
normally leads to passivation. The experimental work has shown that the metallization 
does not involve the surface silicon dimers, which are passivated upon hydrogen or water 
exposure |l| . It was early speculated that the metallization was caused by the incorporation 
of hydrogen in deeper layers, and more specifically in the troughs of the 3x2 reconstruction. 
Recent theoretical works 0, Q, E ] have shown that this is a plausible (although not thoroughly 
accepted) explanation for the surface metallization. Fig. [1] shows the band structure of 
the proposed configuration (schematically shown in the inset), obtained from a Density 
Functional Theory (DFT) calculation [6| in a thick slab containing the top Si layers plus 12 
C-Si bilayers. The shadowed area shows the projection of the bulk band structure, which 
is positioned with respect to the surface bands by aligning the electrostatic potential of the 
bulk with that of the slab calculation at a depth where this potential becomes flat. The 
surface is clearly metallic, with several band crossing the Fermi level. 

As discussed by other authors 0, Q, 5], while some of the metallic bands correspond to 
surface states, others come from states deriving from the conduction band of the host SiC 
crystal. These states can be seen in Fig. [T] as parabolic bands split off the conduction band 
minima at the J and F points of the Surface Brilloiun Zone (SBZ), while the extended con- 
duction band states remain above the Fermi level. In this letter, we study these interesting 
and quite unexplored aspects of this surface, and show how these states become localized 
near the surface forming quantized quasi-two dimensional states due to band bending effects. 
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The key to understand surface metallization is the realization that incorporation of hy- 
drogen to the in-trough sites leads to transfer of electrons from the Si-H-Si center to surface 
and conduction bands that were unoccupied prior to hydrogen incorporation. The separa- 
tion of charges leads to electric fields which produce a strong band bending near the surface. 
Fig. |2](a) (full line) shows the averaged electrostatic potential [7] obtained in our DFT slab 
calculations as a function of the depth. These results show the presence of electric fields at 
depths of more than 10 A from the Si-rich region, well into the host SiC material. This band 
bending acts as a quantum well potential that produces quantization of the conduction band 
states, which become localized in the direction perpendicular to the surface. Fig. [3] shows 
these states, as obtained in the DFT calculations. 

To provide a deeper understanding of the band bending and the appearance of quantized 
states, we present a simple model which captures the physics of the system. We use the 
method proposed by Stern and Howard s|, to study semiconductor inversion layers. We 

only consider the filfirtrnns HnnatfiH tn trip rnnHiiftinn hand hv tfifi in-trmiph adsorbed hy- 
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FIG. 2: (a) Average potential seen by the electrons from the DFT slab calculations (full line) and 
the effective mass model (broken line), (b) z-dependent part of the bound wave functions obtained 
from the effective mass model. Dashed and dot-dashed lines: bound solutions for the longitudinal 
valleys; full line: bound solution for the transverse valley. The origin in the z axis is chosen at the 
position of the Carbon layer closest to the surface. 
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FIG. 3: Constant contour plots of the DFT wave functions of the quasi- two dimensional states at 
the J point (panels (a) and (b)), and the V point (panel (c)). Blue and red indicate positive and 
negative values, respectively. 



drogen atoms. The rest of the electrons (including those donated to surface bands) are not 
considered explicitly, since they do not contribute significantly to band bending effects. We 
work within the effective mass approximation, using the effective mass tensor corresponding 
to the bottom of the conduction band. Since we are dealing with the (001) surface orienta- 
tion, the conduction band has three minima in the SBZ of the lxl surface reconstruction. 
The minima correspond to the conduction band valleys at the X points of the bulk /3-SiC 
Brillouin zone, giving three ellipses with longitudinal and transverse masses of 0.68 and 0.25 



mo, respectively [10]. Two of the ellipses have their longitudinal mass in the surface plane, 
and are centered at the M points of the /3-SiC(001) lxl SBZ. The third one has its longitu- 
dinal mass along the direction perpendicular to the surface, and is centered at the Y point 
of the ^-SiC(OOl) lxl SBZ. Note that the V and M points of the lxl reconstruction fold 
down to the T and J points of the 3x2 reconstruction, which is the reason why the minima 
of the SiC bulk conduction bands projected in the 3x2 SBZ appear, at the same energy, at 
r and J in Fig. □ 

By the use of the effective mass approximation, the problem is reduced to finding the 
solutions of the potential in the direction perpendicular to the surface, while in the parallel 
directions the solutions are plane waves with parabolic dispersion. We model the effective 
potential V e ff(z) seen by the electrons as the contribution of two terms. The first one rep- 
resents the confining potential at the surface, which keeps the electrons inside the bulk. 
Consistently with the effective mass approximation, this potential is taken as zero inside the 
bulk and to raise smoothly at the surface to reach the vacuum level V vac outside the solid. 
Here, V vac is the energy required for an electron at the bottom of the conduction band to 
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escape the surface, which we estimate from the difference between the work function of the 
clean surface and the energy gap of SiC, resulting in a value of 3.5 eV. The second term 
describes the electrostatic attraction between the conduction electrons and the net compen- 
sating positive charge at the surface layers, as well as the repulsion between the conduction 
electrons. The latter is described in the Hartree approximation [8|, thus neglecting exchange 
and correlation effects. The electrostatic term is related to the charge density of the conduc- 
tion band electrons and the associated positive charge left in the adsorbed hydrogen atoms 
by Poisson's equation. 

We solve numerically the ID Schrodinger effective mass equations relating the potential 
V e ff(z) and the wavefunctions Q(z) and their energies for each of the conduction band 
valleys of the SiC(OOl) surface. Here, a = I or t indicates the longitudinal or transversal 
valley, and i is an index labeling each of the bound solutions. The bound states are occupied 
up to the Fermi level, which is determined by the number of electrons donated by the 
adsorbed hydrogen atoms. In all the calculations, the electronic temperature was set to 
zero, so the electronic occupations follow the step function distribution (i.e., there are no 
thermally excited carriers). Once the occupations of the bound states are determined, the 
charge density is constructed and Poisson's equation is solved numerically to obtain the 
updated effective potential V e ff(z). This procedure is repeated until self-consistency is 
achieved. 

The self-consistent electron potential obtained from the model is shown in Fig. [2(a). It is 
very similar to the DFT potential, which clearly shows that the origin of the band bending 
is the donation of electrons to the conduction band. For the quantized wave functions, we 
obtain three bound solutions. For the ellipse perpendicular to the surface, we obtain two 
bound states, at energies 0.34 and 0.01 eV below the bottom of the bulk conduction band, 
shown in dashed and dot-dashed lines in Fig. [2(b), respectively. This ellipse corresponds 
to states at the T point of the 3x2 SBZ. In the DFT results shown in Fig. H} only one 
bound state is observed, with an energy close to that found in the model for the lowest 
state. The higher one is not seen as a bound state in the DFT calculation, a discrepancy 
which is explained by the finite thickness of the slab used in the DFT calculation (note 
that this state extends several tens of Angstrom from the surface, much more than the slab 
thickness). For the ellipses parallel to the surface plane, we obtain only one bound state 
(full line in Fig. E(b)), with an energy 0.16 eV below the bulk conduction band, which is 
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doubly degenerate due to the existence of two equivalent ellipses. It corresponds to states 
in the J point of the 3x2 SBZ. In the DFT calculation, due to the presence of the surface 
reconstruction that breaks the crystal symmetry in the xy plane (due to the presence of 
the dimer rows and troughs), these two states mix up and split. Also, they mix with other 
surface states at similar energies at the J point point. The resulting DFT wave functions 
are those shown in Fig. 0(a)- (b). We note that the lower binding energy and larger spatial 
width of this state compared to the most bound state of the ellipse perpendicular to the 
surface is due to the smaller effective mass of the electrons in the z direction. 

Both in DFT and in the simple model, the underlying bulk SiC crystal is assumed to be 
an intrinsic semiconductor at T=0°K. Therefore, there are no free carriers available from the 
bulk. The band bending effects described here are not due to screening of surface dipoles or 
charges by bulk carriers, as it is the usual case in semiconductor surfaces and interfaces. It is 
purely due to the carriers donated by the adsorbed hydrogen to the conduction bands. For 
substantial doping and finite temperatures, the existence of free carriers will change cuanti- 
tatively the profile of the band bending at long distances, although the qualitative picture of 
band bending originated near the surface from donation of carriers to the conduction band 
remains. For n-type doping, the bulk Fermi level will be also pinned at the bottom of the 
conduction band, so the picture presented here remains accurate. For p-type doping, an 
additional bending will appear at longer distances to align the Fermi level at the surface 
with that of the bulk (pinned around the maximum of the valence band). Its extension will 
depend on the doping level and the temperature, but will be longer than the one obtained 
here, for reasonable values of these variables. 

Finally, it is interesting to note that the computed band bending is able to explain the 



core level shifts experimentally observed by D'Angelo et a/.[ll|. Besides several surface 
related components, these experiments show that the 'bulk' component of the Si 2p core 
level spectra of the clean (unhydrogenated) surface shifts to lower kinetic energies upon 
hydrogenation. The amount of the shift depends on the H exposure, and reaches -0.84 eV 
on saturation. This 'bulk' component of the Si 2p core level spectra is ascribed to the first 
Si layer of the SiC underlying the three top-surface Si planes, because due to the surface 
sensitivity at the photon energy used (150.46 eV) deeper Si layers do not contribute to the 
photoemission signal. The core level shift of this 'bulk' component is therefore a measure 
of the change in the local potential at the position of the topmost silicon layer of the SiC 
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crystal, upon hydrogenation. In our calculations, as shown in Fig. [2j this amounts to -0.93, 
in excellent agreement with the experimental value. 
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